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Abstract 

We analyze lattice data for octet baryon masses from the QCDSF collaboration employing man- 
ifestly covariant Baryon Chiral Perturbation Theory. It is shown that certain combinations of 
low-energy constants can be fixed more accurately than before from this data. We also examine 
the impact of this analysis on the pion-nucleon sigma term, and on the convergence properties of 
baryon mass expansions in the SU(3) symmetry limit. 

PACS numbers: 



I. INTRODUCTION 



Three- flavor Baryon Chiral Perturbation Theory (BChPT) PQ[2] is an effective field theory 
for QCD in the low-energy regime where the lowest-lying meson-octet and the lowest-lying 
baryon octet (and possibly also the decuplet) are the relevant degrees of freedom. 

It is now known for many years that chiral expansions derived from this framework 
are usually not under sufficient theoretical control because of an unsatisfying convergence 
behaviour 1 : higher-order corrections from chiral loops can in general not be guaranteed 
to be much smaller than the leading terms. The quark mass expansion around the three- 
flavor chiral limit converges slowly because the strange quark mass, m s , is much larger 
than the light quark masses mi = m Ui d, and induces sizeable kaon loop corrections to the 
leading tree level results. Detailed discussions on the convergence behaviour of BChPT 
from various perspectives can be found in [3HT0] . The essential point, however, can be 
seen immediately from a comparison of mass scales: while ^2Bq1tli ~ M n w 140 MeV is 
small compared to a typical hadronic scale of ~ 1 GeV, the three-flavor expansions can 
involve y/2B m s =: M Ss « 700 MeV, which is not small compared to that reference scale. 
In fact, loop graphs involving the eta meson are only suppressed by an additional factor of 
Mf)/ (47rFo) ~ 0.5, which can easily be compensated by additional prefactors (Fq is the meson 
decay constant in the three-flavor chiral limit, while Bq is related to the quark condensate 
in the same limit). 

The problem is particularly urgent for chiral extrapolation formulae for lattice simula- 
tions. There, the meson masses are in general even larger than for physical quark masses. 
Consequently, the application of BChPT to such lattice data is not under satisfactory the- 
oretical control. Though one can obtain fits that seem to describe the data fairly well, a 
stable determination of the corresponding low-energy constants (LECs) is not yet within 
reach, because the extracted values can be strongly affected by uncontrolled higher-order 
corrections. For some specific application, one might not even worry much about this, but 
one should recall that the main point of ChPT is to yield relations between many different 
observables, in terms of a limited set (at least at a fixed order in the low-energy expan- 
sion) of parameters. For example, the baryon mass in the chiral limit, mo, and its leading 

1 Of course, we are not concerned here with convergence in the strict mathematical sense, but with the 
reliability of theoretical predictions at low orders in perturbation theory. 
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quark mass corrections will unavoidably appear in practically any complete loop calculation 
in BChPT, in baryon form factors, meson-baryon scattering lengths, etc. Therefore, one 
should try to somehow control the higher-order loop corrections. 

Many alternative ways to overcome this problem have been proposed in the literature, 
e.g. long-distance regularization [6], reordering and resummation schemes 0[TT], large- N c - 
scaling arguments [T2l IT3] and two-flavor expansions for hyperon properties [HJ [15], each 
having its own merits. In this work, we stick to a standard scheme (evaluating the chiral 
loops employing infrared regularization [TU]) and thus analyze the quark mass expansions 
of the baryon octet masses in a fashion closely following [17J. We refer to that reference for 
most of the corresponding technical details (see also [18J). For studies in different schemes, 
see e.g. (H IT9H24"] . We have reproduced the complete one-loop (0(p 4 )) calculation in 
this covariant scheme. Our aim here is to put to use calculations like [17] . by analyzing 
the lattice data presented in [25] which covers a quark mass region where all the masses 
of the pseudoscalar bosons are smaller than the physical eta mass, and to give improved 
estimates for LECs featuring in these calculations. We consider this a first and necessary 
step towards a controlled application of BChPT to such data sets. While in most three- 
flavor lattice simulations, the (large) strange quark mass is fixed 2 , and the physical point 
is approached lowering me (see e.g. [26j EZ] for recent examples), the strategy in [25] is 
different: One simulates at a "symmetric point" where = m s , and then approaches the 
physical point increasing the quark mass difference, but keeping the average quark mass 
fixed. At the symmetric point, the pseudoscalar octet mesons all have the same mass of 
about 400 MeV, and on the trajectory to the physical point, the kaons and etas are lighter 
than observed in nature (there even exist some data points where the meson masses are all 
< 400 MeV at m s = me). This is already a nice feature, but there is a second point that is 
probably even more important. In a surprisingly large region around the symmetric point, 
approximations neglecting all but linear symmetry breaking effects were shown in [22] to 
yield a very good description of the data. In particular, quantities which do not receive linear 
symmetry breaking contributions are observed to be very stable even down to the physical 
point (one example for such a quantity is the famous Gell-Mann-Okubo mass difference). 

2 In [26J, the strange quark mass was 20% larger than the physical value, resulting in large kaon and eta 
masses, which makes the application of three-flavor chiral extrapolations even more problematic. 
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This behaviour can not be understood from ChPT alone, but it puts tight constraints on 
certain combinations of LECs. It is convenient for the purpose of studying the symmetry- 
breaking effects separately to analyze certain baryon mass ratios, where a large factor with 
the dimension of mass (essentially the baryon mass at the symmetric point) cancels out. The 
behaviour of these ratios can be very well described by BChPT (see our Fig. [5j and also [28] 
for a corresponding illustration). Here we want to exploit the advantages just mentioned 
for the determination of the LEC-combinations occurring in octet baryon masses, in the 
hope that this will help e.g. when analyzing lattice data for other quantities with the same 
strategy. 

This paper is organized as follows: After briefly introducing some basic notation in sec. [TX], 



we collect the input from the purely mesonic sector in sec. Ill in some detail. This is 
necessary, because the quark-mass dependence of the masses of the pseudoscalar mesons is 
needed to rewrite the quark mass expansions of the baryon masses at 0(p 4 ) in terms of the 
meson masses. We also include some simple examples of our strategy of analysis in this 



section. In sec. IV, we consider the result of the leading one-loop (0(p 3 )) calculation, to 
set the stage for the later developments. We also show some examples of a first numerical 
analysis, analyzing the experimental baryon masses. Sec. [V] completes the outline of our 
framework specifying the Lagrangians and diagrams entering at fourth chiral order. In 
sec. |VI[ we discuss in detail the extrapolation functions we use, and which combinations of 



LECs can be extracted more accurately than before from the data we consider (and which 
can not). In a series of tables, we present our fit results, obtained with different fit strategies 



and input parameter sets. We discuss our results in sec. |VII| and give a short conclusion of 
our findings. 



II. PRELIMINARIES 



Throughout we work in the isospin limit and set the £ight quark masses to m u = m& =: 
nig. In [25], the quark masses are varied along a trajectory T in the (me, m s )-space where 
the average quark mass m := |(2m^ + m s ) is kept fixed, while 5mg = nig — m is varied. The 
aforesaid trajectory T connects the symmetric point, where rrig = m s = m, and the physical 
point where the strange quark mass attains its physical value while me = \(m u + m^). 
To get familiar with the notation, let us consider the leading terms in the quark-mass 
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expansion of the squares of the pseudo-Goldstone-boson masses. We have 

M 2 := 2B Q me = 2B rh + 2B 5m e , (1) 
M| ;= B Q (me + m s ) = 2B m - B 5m e , (2) 
Ml := ^B (m e + 2m s ) = 2B m - 2B 6m t . (3) 

Here B is related to the quark condensate in the chiral limit, see e.g. [29]. The leading 
term of the pseudo-Goldstone-boson mass at the symmetric point is then found by setting 
Bo5rri£ = 0, 

Ml = 2B Q m . (4) 

Our notation is set up such that quantities in the chiral limit will be denoted by a subscript 
0, while quantities at the symmetric point on T are marked with a * . 
In [25], it was noted that some specific combinations of hadron masses are very stable when 
varying the quark masses along T, compare e.g. Tab. 2 and Fig. 13 in that paper. Below, 
we list three of these combinations which are relevant for our present study: 

X N = \{m N + m 2 + m~) w 1150 MeV, (5) 
Xl = \{2M 2 K + Ml) w (412 MeV) 2 , (6) 
K = \ {Ml + M 2 ) w (400 MeV) 2 . (7) 

The experimental values are given in round brackets. 



III. MESON MASSES AND DECAY CONSTANTS TO ONE LOOP 

Meson masses 

The one-loop expressions for the masses of the pions, kaons and the eta can be found e.g. 
in [29]. We give here the expansions of these formulae in terms of fh and 8mt. The meson 
mass squared at the symmetric point is given (to one-loop accuracy) by 

Ml = 2B rh 1 + - " 128vr 2 (-3L 4 - L 5 + 6L 6 + 2L 8 ) + - log ' 



(4vrF ) 2 V 3 b V A* 
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The expansion of the meson masses in the symmetry breaking parameter 5rri£ reads 

Ml = Ml + (2B 6m e ) (l + Q + 2 log (^^) - 128vr 2 (3L 4 + 2L 5 - 6L 6 - 4L 

'5 + 8 log (^f 2 ^) - 768tt 2 (L 5 - 2L 4 x 



+wm t) ^ — v 24 ; 2Fq2 

(9) 

Mi = M 2 + (2E ^) (-1 - Q + log (^^) - 647T 2 (3L 4 + 2L 5 - 6L 6 - 4L 8 
+ (W ^ J , 

(10) 

Ml = Ml - (2B 5 mi ) (l + Q + 2 log (^^) - 1287T 2 (3L 4 + 2L 5 - 6L 6 - 4L g 

/8 + 12 log f ^f^) + 7687r 2 (L 5 - 12L 7 - 6L 8 ) % 
V w 24tt 2 F 2 



(11) 



Consequently, we have for the combinations X 2 , X 2 : 

J? = \(2Ml + M?) = Aff + <|g^ ( 3 + 41og (^5) - 384^(L 5 - 21.)) , (12) 



(13) 



The low-energy constants Lj of course correspond to the renormalized constants Ll(fi) of 
[29] . In [25], instead of the quark mass difference BoSme, the variable 

, M 2 -X 2 



2 



2Bo5mf f 2Bom 
Ml l 1 + 247r 2 F 2 



+ 



1 + 3 log ( ^ 2B * m ^ j ~ 192tt 2 (3L 4 + 2L 5 - 6L 6 - 4L 8 )^) ^ 

«(7;;:;(^)-i^, i5 - 2i8) ) 



(14) 
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is used to parameterize the SU(3) symmetry breaking. Inverting this relation and again 
expanding up to the second order in the symmetry breaking variable 5rri£, we get 



2B 5 mi = vMl ( 1 - ( l + 3 lQ g ( v/W ° m ) - 192tt 2 (3L 4 + 2L 5 - 6L 6 - AL S 



24tt 2 F 2 V V A* 



(15) 

where we have also truncated the chiral expansion of the coefficient functions in this series. 



Decay constants 



Similarly, the expansion of the decay constants is given by 



'3 + 6 log (^f 2 ^) - 256tt 2 V 



F w = F ic - (B 5m e ) | | + 0((5m e ) 2 ) , (16) 



'3 + 6 log (^f^) - 256tt 2 /V 



F k = F* + (B 5m e ) | | + 0((5m e ) 2 ) , (17) 

'3 + Ql og (VMo^\ _ 2567r 2 L 5 v 
= < ; ■ + ) | | + 0((5m e f) , (18) 



where 

F * = F ° + (iS^p (^ 64 ^( 3L 4 + L 5 ) - 3 log ))+ 0(rh 2 ) . (19) 

Note that if terms nonlinear in the symmetry breaking variable are neglected, the combina- 
tion 

F X :=^(F W + 2F K ) (20) 

is stable at that order. 



7 



Quark masses 



For the sake of completeness, we also give an expression for the quantity M| s := 2B m s 
in terms of the one-loop formulae for the meson masses: 

Ml = 2M 2 K -Ml = 2M 2 K - Ml 

+ ( 3847r2 {M 4 K (±L 4 + 2L 5 - 8L 6 - 4L 8 ) - M*(L 4 + L 5 - 2L 6 - 2L 8 )) ( 21 ) 

+ :>,\i: log (^j - (3M* + 2M 2 M 2 J log ) + 0(M 6 ) . 

Keeping this quantity fixed, one can express the light quark mass dependence of the baryon 
masses through the pion mass dependence, using 

2B m e = Ml (l + U4 l 2F2 (i152tt 2 (M*(2L 4 + L 5 - 4L 6 - 2L 8 ) + M 2 S (L, - 2L 6 )) 

( ^Mj+2Mj^ /^^^ss [ ' 12] 



r2 ,h2 I y^-v^ss j _g M 2 l ( M 



+ (Ml + 2M 2 S ) log 



0(M b 



Numerical analysis 

For an overview on lattice determinations of the parameters Fq, Li we refer to the recent 
discussion in [30]. We reproduce three sets of results in the table below. As central values in 
our numerical analysis, we adopt the set of parameters displayed in the first row (MILC2010). 
The other three parameter sets are used only for the error estimates of our results. The 
numerical values of the Li are given at a renormalization scale of fi = 770 MeV. At the order 

TABLE I: Three parameter sets for the low-energy coefficients Fo,Li. 



F (MeV) 


10 3 L 4 


10 3 L 5 


10 3 L 6 


10 3 L 8 


Ref. 


80.3 ±6.0 


-0.08 ±0.60 


0.98 ±0.40 


-0.02 ±0.40 


0.42 ±0.30 


MILC2010 [31] 


78.3 ±3.2 


0.04 ±0.14 


0.84 ±0.38 


0.07 ±0.11 


0.36 ±0.09 


MILC2009 [32] 


83.8 ±6.4 


-0.06 ±0.10 


1.45 ±0.07 


0.02 ±0.05 


0.62 ±0.04 


PACS2008 [2>] 



we are working, we can replace F in eq. <\8h by F+, see eq. (19). Given the observed stability 



of Att, and the fact that this quantity is even used to set the scale for meson observables in 
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we determine the parameter 2B m from the requirement that M* = X^ hys = 412 MeV. 
We find 



'2B m = M± = (420 ± 40) MeV , (23) 
=> = (112 ±15) MeV. (24) 



We note that inserting the value for F directly in eq. (|8j), one finds very similar values, M+ = 
428 MeV and F± = 112.5 MeV, while the experimental value for F x is Ff ys = 105 MeV. 
As the one-loop effects discussed above enter our calculation of baryon masses formally at 
NNLO (meson masses) or even at two-loop order (decay constants), the estimates given in 



eqs. (23),(|24j) are sufficient for our purposes. 

We also note that the fan plots for the meson masses shown in [25] can be nicely reproduced 
with the values of the LECs employed above, while X n is indeed almost a constant along T. 



For this to be the case, one should expect from eq. (12) that 



/420\ 

3847r 2 (L 5 -2L 8 )^3 + 41ogf — J - 0.58. (25) 

In fact, for our central values (MILC2010) we have 

384tt 2 (L 5 ( / u = 770 MeV) - 2L 8 (/i = 770 MeV)) ~ 0.53 , (26) 

so that terms in X n quadratic in the symmetry breaking are indeed tiny. For X v to be 
equally stable, we infer from eq. (13) that we should have L7 ~ —0.1 ■ 10 -3 . In Fig. [TJ 



we show a "fan plot" for the meson decay constants, i.e. we plot the ratios F n /Fx and 
Fr/Fx as functions of the variable v which parametrizes the SU(3) symmetry breaking, 



see eq. (15). The derivative at the symmetric point (y = 0) is essentially determined by 
the LEC L 5 , which we vary within the error according to MILC2010 to generate the bands 
around the central curves. One observes that the fan plot shows a linear behaviour of the 
ratios with v in the vicinity of the symmetric point, but some curvature is visible near the 
physical point, which seems to be necessary to reproduce the experimental values of F n and 
F K - 



Inserting the experimental pion and kaon masses in eq. (21), one arrives at 



M Ss = (715 ± 100) MeV , (27) 
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FIG. 1: A "fan plot" for the symmetry breaking effects in the meson decay constants Fk and F w . 
The black dots denote the experimental values. 

where again we use the MILC2010 parameters as central values. The uncertainty due to 



the variation of the input parameters is somewhat larger than in eq. (23). We have a clear 
hierarchy 

M n « y/2B m e < y f 2B^n < ^2B Q m s = M- ss , (28) 



compare eqs. (23) and (27). While it is quite obvious that an expansion in M Ss (over a 
typical hadronic scale) will be ineffective, an expansion in M+ can be expected to work 
much better, since the expansion parameter is about 40% smaller. 



IV. BARYON MASSES AT THIRD CHIRAL ORDER 



From the relevant terms in the chiral Lagrangian, 



C+b = {B{iH) - m )B) + ^{B^ l5 {u^B}) + ^(B^ l5 [u^ B}) 

2 2 (29) 

+ b (BB)( X+ ) + b D (B{ X+ , B}) + b F (B[ X+ , B]) + ... 
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one derives the octet baryon masses to 0(p 3 ) in BChPT (see e.g. |19j ) 
m N = m - Ab D M 2 K + 4b F (.M 2 K - M 2 ) - 2b (2M 2 K + M 2 ) 

-6DF to + 2M 3 K - 3M^j + 9F 2 to + 2M\ + mA + D 2 to + 10M| + 9M^) 



96ttF 2 



(30) 



4 . 9F 2 Mi + D 2 to + M| + 3M3 
m A = m + -6 D (-4Mi + M 2 j - 2b [2M 2 K + M 2 ) ±— 2 



(31) 



3 u V ' V u V K ' V 24ttF 2 

D 2 to + 3M£ + Ml) + 3F 2 (M» + 2MI) 
m E = m - 46 D M 2 - 2b (2M 2 K + M 2 J ^ ^pii ^ , 

(32) 

m 3 = m - Ab D M 2 K - Ab F (.M 2 K - M 2 ) - 2b (2M 2 K + M 2 ) 

QDF to + 2M| - 3M^) + 9F 2 to + 2M\ + M z \ + D 2 to 3 + 10M| + 9M 3 ) 

96ttF 2 

(33) 

from the second order quark mass insertions and the sunrise-type loop graph of Fig. [2] To 
this order, the mass of the eta can be eliminated by using the Gell-Mann-Okubo relation, 

3M 2 = 4M 2 - Ml . (34) 

Inserting Mk — > 495 MeV and M n — > 140 MeV gives M v ~ 566 MeV. In our numerical work, 



FIG. 2: Loop graph with axial meson-baryon vertices from the leading-order chiral Lagrangian, 
which yields the leading nonanalytic quark mass correction ~ ^Jm 3 . The full lines represent the 
baryons, while the dashed line represents the pseudoscalar meson tt, K or r/. 

we shall use D — |, F — ^ which is consistent with the available literature [33TI35] . It has 
not been possible so far to obtain these two axial coupling constants from a fit to baryon 
masses. For this choice of D,F, the p 3 -terms cancel completely in the combination 

A DF := - (12m s - 4m A - 5m N - 3m 3 ) = ? {20b D - 3b F ) [m 2 k - M 2 ) + C(p 4 ) . (35) 
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This quantity vanishes in the SU(3) limit, however, it is not too small in the real world. 
Inserting the physical values for the baryon masses (averaging over isospin multiplets), 

m p ^ ys ~ 939 MeV, mf ys ~ 1116 MeV, m p £ ys ~ 1190 MeV, m p 3 hys ~ 1318 MeV, (36) 

gives A P pp S ~ 292 MeV. Assuming that this choice for D, F is roughly correct, and that the 
0{p A ) terms are not too large, one would get y&o — &f ~ 0.65/GeV. It is known, however, 
that the BChPT expressions for the baryon masses only converge slowly at the physical 
point. Moreover, the numerical values of D and F are also a source of uncertainty. We 
shall therefore require only that the estimate for this combination of LECs is valid within a 



— <-b D -b F <-±- forD^^F^ 1 -. (37) 



GeV ~ 3 ~ GeV 4 ' 2 

In contrast to Adf, the well-known Gell-Mann-Okubo difference [36| 137] 

A GM o := | (m E + 3m A - 2m N - 2m H ) = ^ p2 [AM 3 K - M* - 3MJJ + 0(p 4 ) . (38) 

is free of (9(p 2 )-corrections. Employing the input values D — | , F — |, Fq — 80.3 MeV, 
M w = 140 MeV, M K = 495 MeV, M v = 566 MeV yields A GMO = 5.9 MeV + C(p 4 ), while 
Agf/ Q ~ 6.0 MeV with the baryon masses used above. It seems that the higher-order cor- 
rections in this combination are very small: the effect is almost entirely given by the leading 
one-loop correction. The fact that Aqmo is so small numerically could be traced back to 
the fact that it is of order (5mi) 2 in the symmetry breaking variable. 

Fitting the p 3 - expressions to the experimental baryon masses, with D = |, F = ^,rriQ S := 
m — Qb X 2 leads to the results presented in the table below. We vary the input parameter 
Fo appearing in the (9(g 3 )-correction in a wide range to get a first impression of the uncer- 
tainty in these results. The convergence of the expansion of the baryon masses using these 
parameters is poor, however: For example, the third order correction to the nucleon mass is 
—423 MeV for Fo = 80.3 MeV. This is a very large correction, even if one expects a generic 
suppression factor of ~ ~ 0.5 for each additional chiral order. We infer from these 
results that a reliable determination of BChPT low-energy constants at (strange) quark 
masses higher than, or equal to, the value at the physical point is not feasible, due to much 
enhanced non-analytic loop corrections in this mass region. On the other hand, given that 
the above fit results can be regarded as a meaningful first approximation, we should at least 
expect the following pattern to be reproduced also in higher-order calculations: 
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F (MeV) 


mf (GeV) 


b D (GeV- 1 ) 


b F (GeV" 1 ) 


70 


2.279 


-0.039 


-0.912 


80 


2.024 


-0.015 


-0.747 


112 


1.617 


0.024 


-0.484 


140 


1.465 


0.038 


-0.386 



1. lGeV<mf < 3 GeV, 

2. b F is negative, -1 GeV" 1 < b F < GeV" 1 , 

3. bp, is of significantly smaller magnitude than bp. 

Sigma term at third chiral order 

The third-order formula for the pion-nucleon sigma term reads (see also [19]): 

M 2 ■ 4 

a wN {0) = -2M*(26 + b D + bp) - — ^(±o? N M v + 2a%M K + -oP N M v ) . (39) 

where 

a- N = 9 -(D + F)\ a§ = ^5D 2 -6DF + 9F 2 ), a r ' N = \{D - 3F) 2 . (40) 
Numerically, the third order correction gives a contribution of 

at the physical point, to a positive quantity that is of the order 50 MeV ! This again casts 
strong doubts on the applicability of BChPT at the physical point to obtain reasonable 
predictions for sigma terms, without further input. 

V. FOURTH ORDER CONTRIBUTIONS 

At fourth chiral order, additional loop graphs have to be calculated. First, one has a graph 
of the same topology as the graph of Fig. [2j but with an additional quark mass insertion from 
the second order Lagrangian, proportional to b , bp,,F- This graph is shown on the l.h.s. of 
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Fig. [3j In addition, there are tadpole-type contributions with vertices from the second order 
Lagrangian, accompanied by b Q ,bo,F and also eight new constants 6j, i = 1 . . . 4, 8 . . . 11. 
There are also seven new contact terms ~ dj contributing to the masses at fourth order, 
they absorb divergences of the infrared-regularized loop integrals. The self-energy H^ n \p) 
of n-th chiral order is then calculated directly from those graphs, and gives rise to the mass 
shifts Am^ at the different orders: 

Amg } = , Amg } = S (2) , Amf = £ (3) (f> = m ) , 

} d^ 3 \p) ( 42 ) 



=m 



The calculation of the first graph in Fig. [3] is straightforward and does not introduce new 



- . \ u. 



FIG. 3: Loop graphs giving contributions of fourth chiral order. 

free parameters, so we will not discuss it further here (see e.g. Ref. |17j). In the following, we 
are concerned with the contact terms and the tadpole graphs, which involve new unknown 
couplings not present in a simple leading one-loop calculation. These new parameters can be 
grouped in two sets: one which will simply represent new free parameters to be determined 
in the fits, and one which will be fixed using input from other sources. Let us discuss both 
sets in turn. 

Contact terms 

We use the fourth order Lagrangian of Frink and MeiBner [T7] . 

4S = d 1 (B[ X+ , [x+, B}]) + d 2 (B[ X+ , { X+ , B}]) + d 3 (B{ X+ , { X +, B}}) 

(43) 

+ d±(B X+ )(x+B) + d 5 (B[ X+ ,B})(x+) + d 6 (BB)( X+ )(x+) +d r (BB)(xl) . 
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The contact terms give the following contributions to the baryon masses: 
mffi = m ^ t] + B 5m e (A(b D - 3b F ) - 16B m{6d 2 - 44 + 94)) 

- 16{B 5m e ) 2 (94 - 34 + d 3 + 6d 7 ) , 
m W = m [ ct ) + 8B 5mt (b D + 16B md 3 ) - 96(B 5m e ) 2 (24 + 4 + d 7 ) 



(ct) 

= m 



(ct) 

m- = m 



{ ct) - 8B 5m e (b D + 16B md 3 ) - 16(B Q 5m e ) 2 (44 + 64) 
[ c4) + B 5m e (4(b D + 3b F ) + 16B m(Qd 2 + 4d 3 + 94)) 



(44) 



- 16(B Sm £ ) 2 (94 + 34 + d 3 + 64) , 
where m[ ct) = -4B m (3b + 2b D ) - 16(B m) 2 (44 + 9d e + 3d 7 ) . 

We also recall that the combination has no linear symmetry-breaking terms. In partic- 
ular, 

Xf ] = m{ ct) - 32(B 5m e ) 2 (34 + d 3 + 34) . (45) 

In the contact term contributions linear in the symmetry breaking, only the combinations 
bo '■= bo + lGBofnd 3 and bp '■— b F + 4:Bofn(2d 2 + 34) enter. While 4,3 also enter the terms 
quadratic in 5mi, we can not determine 4 from the fan plots at the order we are working 
(the tadpole contribution involving bp is already of fourth order in the chiral counting). 
Moreover, the LECs 60 an d 4 only enter in the combination 

m ^ = m - l2B m (b + 12B md 6 ) + ... , (46) 

so these two parameters can also not be determined individually from the fan plots alone 
(again, the shift proportional to 4 would be of higher chiral order in the tadpole graph 
involving 60) ■ For a later application, we also give the contact term contribution to the 
pion-nucleon sigma term. It is computed here employing the Hellmann-Feynman-theorem, 



2<9mi rf) 1 (lct) 2 ( dm% fit) {2 , n 
= mi- — + -B m e m K N ' + (B dmi)- I mi— — h B Q mern> N ' 



(47) 



Inserting the expressions from eq. (44), we find 

(32 8 
4& + 2b D + 2b F + —B m (44 + 94 + 34) + ^B m (64 - 44 + 94 



+ 16(^0^) (34 + 4 - 4 + 34 + 2d 



7/ 



(48) 
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Tadpole graphs 



The second order Lagrangian we employ here differs slightly from the one given in [T7] 
because since then one further term could be eliminated [38J . The Lagrangian is then given 
by 

C% = b D/F (B[ X+ ,B] ± ) + b (BB)( X+ ) 

+ b 1/2 (B [u^ K, B] T }) + h{B{u„ K, B}}) + b A (BB)(uX) 
+ i(b 5/6 (Ba^ [K, u v ] , B] T ) + b 7 {Ba^ Uil ){u v B)) 

+ ^ ((BY K, K [D u , B}] T ] ) + (BY [A,, K, K, B]] T ] )) (49) 

+ ^ ((BY {u„ K, [D u , B]}}) + (BY [A,, {«", 5}}])) 

+ ^(2(By [D^SJXttX) + (5 7 ^}([ J D i ,,^] U i/ + « /1 [ J D, ) ^]}) . 

This form has also been used in HQ], and we shall also take the numerical values of the 
LECs &1-4, 6s— n from the latter references. Note, however, that in those references, the 
chiral potential derived from the above Lagrangian is iterated to infinite order in a coupled- 
channel Bethe-Salpeter equation, introducing a model-dependent uncertainty not controlled 
within strict ChPT. To estimate the impact of this uncertainty on our results, we perform 
all our fits with three different parameter sets from meson-baryon scattering (set 1-3). In 
addition, we also add one parameter set where we take all &i-n to vanish (set 4). The four 
parameter sets are collected in the table below (in GeV -1 ). 65-7 do not contribute to the 

TABLE II: Four parameter sets for the low-energy coefficients 61-4, &8-11- 



set 


h 


b 2 


bs 


64 


^8 


kg 


bw 


611 


Ref. 


1 


-0.082 


-0.118 


-1.890 


-0.215 


0.609 


-0.633 


1.920 


-0.919 


m 


2 


-0.126 


-0.139 


-2.227 


-0.288 


0.610 


-0.677 


2.027 


-0.847 


|40| 


3 


-0.014 


-0.207 


-1.063 


-1.312 


0.272 


-0.483 


1.054 


0.328 


|40| 


4 





























baryon masses at one-loop order. We shall not write out the full tadpole contributions to 
the baryon masses. The corresponding contributions to the terms linear in the symmetry 

16 



breaking 5rri£ and to the baryon mass in the symmetry limit, m*, can be read off from 
eqs. (69) and (71 ). 



Renormalization 

The infrared regularized loop integrals still contain UV divergences as d — > 4, proportional 

to 

i: =£(dV^) + r ( i ) + i]). Po) 

The couplings &o,d,f do not have to absorb divergences from the infrared regularized loop 
graphs. At third order, there are no counterterms contributing to the baryon masses, and 
therefore also no divergences. The divergences of the infrared singular parts of the loop 
integrals at fourth order can be absorbed by the following renormalization prescription for 
the counterterms df. 

rf l = 7 ! R ^ + < (r) (/i), (51) 

with 



1 / D 2 — 3F 2 

^ = ^2 M14 + Q9D 2 + 81F 2 ) + 162b F DF + (126 x - 46 3 + 36 8 - b w} i , 

(Zt \ m 



72 R = 
12 48 



1 / 6DF \ 
-2 120b D DF + b F (4 + 60D 2 + 108F 2 ) + + 3(46 2 + b 9 ) 



1 / 3( D 2 — 3F 2 ) 

7s R = 66 D (4 + 13D 2 + 9F 2 ) + 108b F DF + -± - (366a + 46 3 + 9b 8 + b 

4o-^o \ m 



2 I uf£)^ -r ioiy -r »-r ; -r luoupiyr -| ^ou^i -r ^3 -r yu 8 t uio 

IR ' 



1 / 9(D 2 -3F 2 ) 
7f = ^2 -4M11 + - J ^ + 108^ - 46 3 + 276 8 - b 10 , , 



72F 2 \ " K ' m 

- (1086a + 140&3 + 2646 4 + 276 8 + 356i + 66611; 



1 / 1 7n 2 _i_ q p2 

7 i R = 1206 + 6 D (28 - 24(7D 2 + 9F 2 )) - 4326 F DF 

144^ q \ m 

- (366i + 6863 + I2O64 + 96 8 + 17&io + 306n) 
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(52) 



We can then define scale-independent quantities di by 

167T 2 ^ fi 



4%)^ + ^^). (53) 



VI. ANALYSIS OF FAN PLOT DATA 

The main objective of this study is the analysis of the fan plot data collected in Tab. 22 
and Fig. 20 of [25]. The following ratios of baryon masses are considered, 

f B :=^ for B = N,A,E,Z, (54) 

see also eq. ^ for the definition of X^. At the point on the symmetric line where mi = 
m s = rh, we have fg = — = 1 and X v = M± = 412 MeV. The symmetry breaking is then 
switched on by varying Srri£ for fixed rh. For the data points with the greatest distance 



from the symmetric point, we compute v = —0.692 from eq. (14). The meson masses at 
this point are M n « 229 MeV, M K « 477 MeV and (using eq. Q) M v « 535 MeV. All 
meson masses stay well below the experimental eta mass, in contrast to other lattice data 
sets (with m s fixed instead of rh), where the meson masses are usually larger than the 
meson masses in the real world. There is therefore good reason to believe that the chiral 
extrapolation formulae provided by three-flavor BChPT will work better when applied to 
the present data set than in the usual applications, though meson masses of ~ 400 MeV are 
probably still too high to assure a proper convergence behaviour of the chiral expansions 
in general. Another indication of a better behaviour of the quantities considered here is 
the observation (made in [2S]) that finite- volume effects, which are mainly generated by the 
chiral loops, tend to cancel out in the baryon mass ratios In this section, we will see 
that certain combinations of low-energy constants can already be extracted quite reliably 
from the fan plot data for the ratios Jb, namely, the LEC combinations which parameterize 
the leading symmetry breaking contributions to the baryon masses. We consider this as a 
first and necessary step towards a theoretically controlled application of BChPT formulae 
to lattice data. We also make an attempt to determine the remaining LECs occuring in the 
baryon mass extrapolations, however, the constants parameterizing the singlet contributions 
can only be roughly estimated. 

Though we shall use the full one-loop BChPT expressions for the baryon masses in the 
ratios fs in the end, it is instructive to consider the expansion of the latter ratios in the 
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variable 5m£ parameterizing the symmetry breaking. This will make clear which subsets 
and combinations of LECs can be more accurately determined than before, and what the 
qualifications and sources of uncertainties for the various determinations are. 



Expansion of mass ratios /# 

Employing the full one-loop ChPT calculation of the baryon masses, we are sensitive only 
up to terms (5me) k fh 2 ~ k , while higher order terms will be modified by terms on the two-loop 
level. Therefore, the following equations are to be understood as resulting from a double 
expansion: the chiral expansion, on the one hand, and the expansion in 8me, on the other 
hand. One finds 

f B = 1 + (B Sm e )f^ + (B 5m e ) 2 f^ + 0((5m e f) , (55) 

where 



4(6d - 3b F ) AT 7T-3( J D 2 + 10DF - 3F 2 ) 

= \/ ZB Q m — -2 

mo 62m irl < Q 

V m o m J 

Hbr, , z(d 2 -zf 2 ) 



,(D _ *b D /^7- 3(D 2 - 3F 2 ) 



and 



+ (1650m 



f (D_ 86 C „/^ 3 (- p2 ~ 3F2 ) 
J* ~ — — + s/lBtfn 2 



m 16m 7iFQ 

- (MHrfft) f 2 ( 3 "° + f°>"° + ?* + . . ) + £>(*«/») , 

V m o m / 

4(6 D + 3b F ) A 77T-3( J D 2 - - 3F 2 ) 

JS = V 2B m ™ ^2 



mg m 



^ 128v^mm 7rF 2 m "' V ; ' 



(56) 



(57) 



(58) 



(59) 



(2) _ D 2 + 18DF-3F 2 3^-34-4 1/2 
/jV " 128v^mm 7rF 2 16 m + ■ ■ ■ + U(m ) , 

/f = - "(g^) + 32 3 ^- 5 ^- 3 ^ + . . . + 0(mV 2 ) , (61) 
128v / 25o^m 7rF 2 m 



(62) 
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128y/2B mm 7iFQ m 
The dots stand for corrections from loop graphs of fourth chiral order, which we do not 

display explicitly here. Let us concentrate first on the contribution from the contact terms. 

We make the following observations: 

i f (i) _ f (i) 

2- f$> + + fi n) = 0, 

3. the baryon mass m only appears in combination with low-energy constants, in the 
ratios D/^/mo, F/^/mo, bi/m , di/m , 

4. the LEC d 7 does not appear at this order in the double expansion of the ratios. 

The last two items would lead to problems when only using the fan plots and ratios fs in 
the fitting procedure: First, since terms of order (^m^) 3 and higher are expected to be small 
in most parts of the fan plots, it is very likely that a stable determination of d-j would be 
difficult. Second, we would like to use the values of D, F, &i_4, b$-u directly as input, and 
the baryon mass in the chiral limit, mo, is not accurately known. What is more or less 
accurately known, however, is the baryon mass at the symmetric point, 

m* = m - 4^0771(360 + 2b D ) + C(m 3/2 ) , (64) 

for reasons discussed in the next subsection. We can eliminate mo in favor of m*. Moreover, 



we can replace Fq by (see eq. (19)) at the order we are working. These replacements are 
rather natural when expanding in the symmetry breaking around the point me = m s = m. 
Doing this, the previous expansions read after this slight reordering (in terms of bn,F defined 
above) : 



/£> = ^ UD - ° UF> - 



4(6d - 3b F ) /^-3( J D 2 + 10DF - 3F 



2\ 



m- 



m* 32m i ,irF 2 
AB m \ b D -3b F /5 | sb ( V2B^rh \\ 
\3 \ ji J ) 



(4vrF, 

+ ^(13D 2 - 30DF + 9F 2 ) (-+log( ) (60) 

m* \3 \ fi ' ! 



(5D 2 - QDF + 9F 2 ) - + log . 
m* V 3 V A* 

D 2 + 10DF-3F 2 (3 , fVW^rh 
+ a [-, + log — 
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+ + (9(m 3 / 2 ) 



f(l) 

J a 



8b 



D 



'2B m 



-3(D 2 - 3F 2 



IQm+irF? 



46 



AB m 
~ (4ttF,) 2 

V 



26 



m 4 



+ 



5 
3 

726 F 



+ 8 log 
DF 



V A* 



+ log 



(66) 



+ 



2(D 2 - 3F 2 



mf 



+ log 



T A (1) + 0(m^) 



f(l) 



J A ' 



(1) 

A. 

(h _ / ,(1) 



M' = -(f ( N> + 



(67) 
(68) 



while the formulae for f B just change in that m — > m* at the chiral order we are working. 

stand for contributions from tadpole graphs proportional to 61-4, &s-ii- These are 
given below. We see that from the derivative of the ratios at the symmetric point, we can 
determine two parameters only, due to the symmetry constraints in (i) and (ii), namely &d,f- 
From the tadpole contributions to the 



T 

1 N 



(1) 



T 

'A 



(1) 



4B m 
3(47rF,) 2 

8B m 
3(4vrF,)2 



9&i - 15& 2 + 63 



96i 



&3 



1 + 4 log 



+ 



96« - 156q 



'10 



1 + 4 log 



/J 



+ 



96« 



no 



log 



log 



(69) 



we can see that 62,9 enter there only for N, S, while the combinations 9&i + 63, 9&8 + &10 enter 
for all four baryons. 



Expansion of Xn 

The above ratios fs are all normalized to Xjy := \{jn^ + mj + 777.5). This combination 
was shown in [25] to be practically constant along a trajectory with constant rh, and thus 
approximately equal to m™ um = 1150 MeV for the choice of 777 corresponding to the experi- 
mental hadron masses. Subsequently, was even used to set the scale in the simulations 
leading to the fan plot for the baryon masses. Employing our ChPT formulae, the flat 
behaviour of Xn is not automatically guaranteed, but should be enforced in the fitting pro- 
cedure. To see how this can be done, we write down the expansion of Xm in the symmetry 
breaking variable Smf 

X N = m+ + (B 5m e ) 2 xff + 0((5m e ) 3 ) , (70) 
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with 



f5£) 2 + QF 2 ) 

m* = m - AB m(3b + 2b D ) - (2B mf 2 K ^-^ ; 



^2 



+ 



(2B m 
3(4ttF,) 2 

- 2(96 8 + 



32(36 + 26 D )log 



V A* 



8(96! + 76 3 + 126 4 ) log 



76io + 126 u ) (log " ^ - 12(4ttF + ) 2 (4rf 3 + 9d 6 + 3rf 7 ) 



(71) 



_(5D 2 + 9F 2 ) / i + 2lQ /yggg\ 
m V \ V J 



+ G(m 5 / 2 ) , 



and 



(2) 



A? 



3(D 2 + 2F 2 ) 



- 32 (3di + d 3 + 3d 7 



5b 

> F 2 



7T 



+ 



2D 2 5 + 6 log 



6vr 2 F 2 
(126! + 8fe 3 + 156 4 ) 

12vr 2 F 2 
(126 8 + 8610 + 156n 

24tt 2 F 2 
D 2 + 2F 2 /5 



U/2B^n 



15 + 20 log 



. N2B^n 

log 

V V 



3 + 41og 
l + 21og 



f y/2B^rh 



f ^2B^ rn 



(72) 



2m 7r 2 F 2 \4 



+ 



rn 



l/2\ 



Here, the LEC d 7 shows up at order m 2 and (5m^) 2 , respectively. The baryon mass in 
the chiral limit, m , and the combinations of LECs appearing here will be determined in 
the next section from the behaviour of m*(m), for which some lattice data points with 
300 MeV < M* < 500 MeV exist. We repeat that we dismiss any lattice data involving 
meson masses much above 500 MeV. 



Linear fits to the fan plots 

As a first step, we perform a linearized fit to the fan plot data only. Only the LEC- 
combinations that enter terms linear in the symmetry breaking can be safely determined 
from such fits. As shown above, in a full one-loop calculation in ChPT, these are the 
combinations &d,f- As an additional constraint, we impose that the mass combination Xn 
does not deviate much from the physical value ~ 1150 MeV. Approximating 

w + ,/£-" . m 

X N ^m* + v 2 xf u) , 
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using eqs. (68) and (15), and truncating the coefficients after O((2B m) 2 ), we can determine 



the following four combinations of LECs from such fits: 

&d = bp + 16B Q rhd 3 , bp = bp + 4B m(2d2 + 3d 5 ) = b' F + 8B md2 , 

1 (74) 
— ( 
3 

The loop functions are evaluated at a scale fi = 1150 MeV. We have checked that, given the 



b = b - 2B m ( -d 3 - 6d 6 - 2d 7 j , d 7 = di + -d 3 + d 7 . 



running of the coupling constants with /z from eq. (53), the masses are independent of the 
choice of scale when truncating after 0(p 4 ). 

In our numerical analysis, we use 2B^m = (420 MeV) 2 in order to fix the meson mass at 
the symmetric point, M*, to 412 MeV for the set of meson LECs from Tab. [TJ m* is set 
to rn™ um = 1150 MeV in the expressions for fg' u \ eliminating the dependence of these 



quantities on mo (compare eq. (68)), while F* is varied in the range (112 ± 30) MeV to get 
an estimate for the error of the determined combinations (one finds that the results are not 
very sensitive to the choice of 2Bofn). The results of these fits are given in the tables below, 
for some given values of mo and F*, and for all four parameter sets of &i-n from Tab. [XTJ 
Based on the experience with these fits, and on the expectation of higher order corrections to 



the approximations made in eqs. (73 ), we can impose the following bounds on the parameter 



combinations in question (in appropriate units, specified in the above tables): 

0.05 <b D < 0.15, (75) 

-0.50 <b F < -0.25, (76) 

-1.50 <b < -0.20, (77) 

-0.35 <d 7 < -0.05. (78) 

The parameter combinations bp,,F will be fixed in the full fits for every set 1 — 4, while 
mo, bo, d\- 7 will be left free - the comparison with the values bo, d 7 from the tables above will 
only serve as a consistency check afterwards. Please note that the numerical values given in 
the previous four equations relate to the fixed values for \i and 2B m specified above. 



Full one-loop fits to the fan plots 



In the next step, we use the full one-loop expressions for the baryon masses in the 
fit functions. We fit to the fan plot data for the baryon octet masses (Table 22 of 
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TABLE III: bi-u from set 1 



F* (MeV) 


m (MeV) 


b D (GeV 1 ) 


b F (GeV" 1 ) 


bo (GeV" 1 ) 


d 7 (GeV" 3 ) 


112 


800 


0.078 


-0.352 


-0.888 


-0.158 


112 


1000 


0.078 


-0.352 


-0.764 


-0.168 


112 


1200 


0.078 


-0.352 


-0.639 


-0.180 


80 


800 


0.109 


-0.434 


-1.158 


-0.257 


80 


1000 


0.109 


-0.434 


-1.067 


-0.271 


80 


1200 


0.109 


-0.434 


-0.974 


-0.287 


140 


1000 


0.066 


-0.306 


-0.599 


-0.118 



TABLE IV: &i_n from set 2 



(MeV) 


m (MeV) 


b D (GeV" 1 ) 


b F (GeV" 1 ) 


bo (GeV- 1 ) 


d 7 (GeV" 3 ) 


112 


800 


0.062 


-0.355 


-0.961 


-0.173 


112 


1000 


0.062 


-0.355 


-0.838 


-0.183 


112 


1200 


0.062 


-0.355 


-0.713 


-0.195 


80 


800 


0.088 


-0.437 


-1.270 


-0.277 


80 


1000 


0.088 


-0.437 


-1.179 


-0.291 


80 


1200 


0.088 


-0.437 


-1.086 


-0.307 


140 


1000 


0.054 


-0.308 


-0.652 


-0.130 



[25] ) and the lowest three points for Xj^(m,u = 0) = m*(M*) (from Table 19 of 
[25] ) where M* < 412 MeV. In addition, we require that X N (v) « 1150 MeV for 
v e {-0.692,-0.558,-0.404,-0.275,-0.128,0,0.181} (i.e. at the fan plot data points), 
where we allow for an error of 10% for X^{v) at those points. m*(M*) is derived from 



eq. (71) with the help of eq. (|8j). Throughout we eliminate F in favor of F+ - the differ- 
ence, however, formally amounts to an 0(p 5 )-eSect in the baryon masses, which is beyond 
the order we are working in. The ratios fs are directly obtained by inserting the one-loop 
expressions for the baryon masses, without further expansion in v or rh. Also, Xn{v) is 
not truncated after Oiy 2 ) in the full fits. In the fan plot fitting functions, including Xn(v), 
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TABLE V: &i_n from set 3 



F* (MeV) 


m (MeV) 


b D (GeV" 1 ) 


b F (GeV" 1 ) 


bo (GeV- 1 ) 


d 7 (GeV~ 3 ) 


112 


800 


0.083 


-0.358 


-1.024 


-0.181 


112 


1000 


0.083 


-0.358 


-0.900 


-0.192 


112 


1200 


0.083 


-0.358 


-0.775 


-0.203 


80 


800 


0.117 


-0.443 


-1.362 


-0.289 


80 


1000 


0.117 


-0.443 


-1.271 


-0.303 


80 


1200 


0.117 


-0.443 


-1.178 


-0.319 


140 


1000 


0.069 


-0.310 


-0.697 


-0.136 



TABLE VI: &i_n from set 4 



(MeV) 


m (MeV) 


b D (GeV" 1 ) 


b F (GeV" 1 ) 


bo (GeV- 1 ) 


d 7 (GeV" 3 ) 


112 


800 


0.072 


-0.305 


-0.504 


-0.091 


112 


1000 


0.072 


-0.305 


-0.380 


-0.101 


112 


1200 


0.072 


-0.305 


-0.255 


-0.113 


80 


800 


0.095 


-0.365 


-0.580 


-0.166 


80 


1000 


0.095 


-0.365 


-0.490 


-0.180 


80 


1200 


0.095 


-0.365 


-0.396 


-0.196 


140 


1000 


0.063 


-0.271 


-0.322 


-0.068 



where the average quark mass is fixed to its physical value, we shall use the numerical val- 
ues from eqs. (23 24). The meson-LECs will be taken from the MILC2010 set, see Tab. |TJ 
The two combinations bp>,F will then be fixed to the corresponding values determined in the 
previous subsection, for each of the four sets of &i_n from Tab. [TTl 



With given values for bp>,F for y/^Bofn = 0.42 GeV, we can determine seven more param- 
eters from the present data set. First, the (<5m£) 2 -terms in the baryon mass formulae can be 
fixed by four free parameters, which may be taken as 

1 4 

d\ := d\ - -d 3 , d 2 , d^:=d4 + -d 3 (79) 
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see e.g. eq. (63)) and d 7 denned in eq. (74). Furthermore, inserting 



b D = b D - 8(0.42 GeV) 2 d 3 (80) 



in eq. ( 71 ) , we can determine the combinations 
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&;:=&o--(0.42GeV) 2 ci3, 

-41 < 81 > 
d 6 := -d 3 + d 6 + -d 7 

and rriQ from the running of X^iSTni = 0) with B^m. The combination bo used in the 



previous subsection can then be computed, for y/2B m = 0.42 GeV, as 

bo = b' + UBofhde . (82) 



The results are given below, in Tab. VII The fits are very good, % /d.o.f. ~ 0.2. This can be 
understood from the fact that the fan plots look very linear near the symmetric point, and 
the linear approximation fixed by the parameter combinations bo,F already describes this 
part quite well - better than one could have expected beforehand. Moreover, we have enough 
free LECs at hand to describe the behaviour of in a satisfying manner. The emerging 
fitting parameters are all of natural size, and itlq ~ 1 GeV. In the fits leading to the results 



in Tab. VIII, we have also included the physical (experimental) values for the baryon masses 
in the fit, assigning an error of 5MeV to each due to our neglection of isospin-breaking 
effects. The baryon masses at the physical point are obtained inserting v = —0.885 and our 
standard value for 2Bofh in our mass formulae. As a consistency check of our procedure, we 
have also performed the fits without the constraints on bo.F- The corresponding results can 
be found in Tab. |IX| and Tab. |X| To get an estimate of the stability of our results due to 
uncertainties in the input values and higher order corrections, we have also performed full 
fits where the meson decay constant is strictly truncated to the value F of the MILC2010 
set of LECs, see Tab. |T} The fit parameters following with this input can be read off from 
Tab. |XI| We observe a good overall agreement of all these fits: varying the input parameters 
bi, F+, Li always leads to similar results. We will discuss our results in detail in the next 
section. 
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TABLE VII: Results for the fit parameters. The set-number refers to Tab. |ll[ bp> t p are fixed input, 
mo is given in GeV, b ,}>D,F in GeV -1 , and the di are given in GeV -3 at a scale 11 = 1150 MeV. 
Here, the experimental masses are not included in the fit. 



set 


b D 


b F 


m Q 


b'o 


dt 


d 2 


C?4 


d e 


d 7 


bo 


1 


0.078 


-0.352 


1.029 


-0.364 


0.035 


0.122 


0.020 


-0.532 


-0.244 


-0.927 


2 


0.062 


-0.355 


1.051 


-0.373 


0.040 


0.126 


0.018 


-0.603 


-0.272 


-1.011 


3 


0.083 


-0.358 


1.064 


-0.400 


0.030 


0.121 


0.030 


-0.641 


-0.289 


-1.078 


4 


0.072 


-0.305 


0.930 


-0.271 


0.025 


0.101 


0.042 


-0.213 


-0.113 


-0.497 



TABLE VIII: Results for the fit parameters. Here, the fit includes the experimental baryon masses. 



set 


b D 


bp 


m 


K 


di 


d 2 


C?4 


4 


d 7 


bo 


1 


0.078 


-0.352 


1.026 


-0.370 


0.035 


0.046 


0.044 


-0.525 


-0.257 


-0.926 


2 


0.062 


-0.355 


1.046 


-0.381 


0.041 


0.050 


0.041 


-0.594 


-0.287 


-1.010 


3 


0.083 


-0.358 


1.058 


-0.411 


0.030 


0.047 


0.054 


-0.630 


-0.306 


-1.077 


4 


0.072 


-0.305 


0.936 


-0.262 


0.024 


0.017 


0.063 


-0.222 


-0.115 


-0.496 



VII. DISCUSSION 



Convergence at the symmetric point 



We start the discussion of our results with the singlet sector where 5me = 0. Let us have 



a look at the expansion of m*, using the parameter sets from Tab. VII as typical examples. 
The dimensionless suppression factor could naively be expected to be of the order 

y/2i%rh 420 



4ttF* 



4tt • 112 



0.3. 



(83) 



Of course, we already know, e.g. from the considerations in sec. IV, that the leading nonan- 



alytic loop correction ~ M 3 is much enhanced w.r.t. this expectation. Using eq. (71), we 
find 



m - 4B m(3b + 2b D ) - (2B m 



(5L> 2 + 9F 2 ) t {2B rh) 2 
' 3(4ttF*) 2 



24ttF= 



+ 



+ 0(m 5 / 2 ) 
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TABLE IX: Consistency check. Here, the combinations bp>,F are not fixed, and the experimental 
masses are not included in the fit. The values for bo, bp are computed for each set of fit results. 



set 


m 


b'o 


b D 


b' F 


di 


d 2 


di 


d 6 


d 7 


k 


b F 


1 


1.028 


-0.372 


0.087 


-0.405 


0.043 


0.119 


0.039 


-0.531 


-0.244 


-0.933 


-0.321 


2 


1.050 


-0.381 


0.072 


-0.409 


0.048 


0.122 


0.036 


-0.602 


-0.273 


-1.017 


-0.323 


3 


1.063 


-0.407 


0.092 


-0.412 


0.040 


0.121 


0.051 


-0.640 


-0.289 


-1.084 


-0.327 


4 


0.929 


-0.278 


0.079 


-0.346 


0.033 


0.095 


0.059 


-0.211 


-0.114 


-0.501 


-0.280 



TABLE X: Consistency check. Here, the fit includes the experimental baryon masses, but the 
combinations br>,F ar e not fixed. 



set 


mo 


^0 


b D 




di 


d 2 




d 6 


d 7 


^0 


b F 


1 


1.021 


-0.378 


0.077 


-0.357 


0.039 


0.022 


0.051 


-0.518 


-0.257 


-0.926 


-0.341 


2 


1.042 


-0.389 


0.062 


-0.361 


0.044 


0.025 


0.048 


-0.587 


-0.287 


-1.010 


-0.343 


3 


1.054 


-0.417 


0.082 


-0.365 


0.035 


0.025 


0.062 


-0.624 


-0.305 


-1.077 


-0.347 


4 


0.930 


-0.269 


0.069 


-0.295 


0.027 


0.010 


0.067 


-0.213 


-0.115 


-0.494 


-0.302 



(1029 + 330 - 397 + 186) MeV = 1148 MeV for set 1 from Tab. |VH 



(1051 + 351 - 397 + 143) MeV = 1148 MeV for set 2 from Tab. EH 
->• (1064 + 365 - 397 + 115) MeV = 1147 MeV for set 3 from Tab. |VH 
-> ( 930 + 236 - 397 + 380) MeV = 1149 MeV for set 4 from Tab. IVTT1 

The result for m* is very close to rn™ m in all cases due to the constraint included in the 
fit. The coefficient in square brackets giving the fourth order contribution can be read 



off from eq. (71). We see that the observed convergence behaviour is still inconclusive, 
though the fourth order contributions employing the LECs &i_n from the meson-baryon 
scattering amplitudes are roughly in line with a suppression factor of ~ 0.5. A similar 
pattern (at physical quark masses) has e.g. been found in [21], where some fourth order 
effects were estimated. In contrast, the fits for &i_n = result in expansions of which 
seem meaningless in the sense of perturbation theory. 
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TABLE XI: Stability check. Here we use the value from MILC2010 for the meson decay constant 
as input everywhere, instead of F±. The fit includes the experimental baryon masses, but the 
combinations bp> p are not fixed. 



set 


m 


b'o 


b D 


b'p 


di 


d 2 


di 


d 6 


d 7 


^0 


bp 


1 


1.116 


-0.329 


0.118 


-0.470 


0.075 


0.121 


0.134 


-1.247 


-0.524 


-1.649 


-0.385 


2 


1.144 


-0.337 


0.097 


-0.476 


0.084 


0.129 


0.123 


-1.411 


-0.583 


-1.831 


-0.385 


3 


1.161 


-0.371 


0.126 


-0.482 


0.067 


0.127 


0.157 


-1.500 


-0.618 


-1.960 


-0.392 


4 


0.993 


-0.211 


0.098 


-0.381 


0.048 


0.063 


0.152 


-0.516 


-0.250 


-0.757 


-0.336 




0.2 0.3 0.4 

M* [GeV] 



0.6 



FIG. 4: The function m+(M+) for the four parameter sets from Tab. VII| (full lines) and from 



Tab. XI (dashed lines). Red: set 1, green: set 2, blue: set 3, black: set 4. The band shown here is 
the one-sigma error band pertaining to the full red curve. 

To further study the convergence issue, we set the average quark mass to 2Bom = 



(300 MeV) 2 (this entails 104 MeV using eq. (19)). Then we would have 



->- (1029 + 169 - 168 - 10) MeV = 1020 MeV for set 1 from Tab. |VH 
->- (1051 + 179 - 168 - 39) MeV = 1023 MeV for set 2 from Tab. 
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m* -> (1064 + 186 - 168 - 57) MeV = 1025 MeV for set 3 from Tab. EU 
m* (930 + 120 - 168 + 122) MeV = 1004 MeV for set 4 from Tab. IVTil 



The third order term is still comparatively large, cancelling most of the second order 
corrections, but the fourth order contributions already seem to be well under control for 
this low average quark mass. Unfortunately, there is no data available for such low average 
quark masses yet. 

The parameters which are mostly determined from the behaviour of the singlet sector 
can therefore not be fixed very accurately from the data set considered here. In particular, 
mo, Vq and d§ will be subject to large uncertainties. This also reflects itself in the following 
experiment: Including also the fourth data point for m*(M+) at M± ~ 530 MeV in the fit, 
we observe that the aforesaid three parameters vary rather drastically (e.g. the results of 



Tab. |VIl| shift to m ~ 800 MeV, b' ~ -0.9 GeV -1 , d e ~ 0.04 GeV" 3 ), while the remaining 
parameters remain roughly the same. This is particularly disturbing for mo, a constant that 
reappears everywhere in three-flavor chiral extrapolations for baryon observables. We would 
like to point out that it would be very helpful for this purpose to have some more m^-data 
points for M+ < 300 MeV. Up to now, we can only give the following very rough bounds, 



based on Tab. VII XI and the variation of the input parameters: 



800 MeV <m < 1200 MeV , (84) 
-1 GeV" 1 <b' < GeV" 1 , (85) 
-2 GeV" 3 <d 6 < 0.50 GeV" 3 . (86) 



As a consequence of the results in Tab. XI, we revise the lower limit for the allowed range 



of the combination 6 (see eq. (77)) 



-2 GeV" 1 <b < -0.20 GeV" 1 . (87) 

In Fig. |4| we show m* as a function of M*. One immediately sees that the uncertainty grows 
rapidly above the region where the curves are fixed by the fit to data (note that the fourth 
data point shown here is not included in the fits). Recall that the only difference between 
the fits leading to the full and the dashed curves in this figure is the C(p 5 )-effect F± — > F . 
Beyond M* > 500 MeV, this higher-order effect has a very large impact on m*. Of course one 
could obtain good fits to lattice data in this region with both choices for the decay constant, 
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but in our opinion, these fits would not be reliable in the sense of a stable determination of 
the fitting parameters (LECs) . The finding that chiral extrapolations should not be trusted 
for meson masses much above 500 MeV is consistent with related studies [HI SI]- Beyond this 
regime, the extrapolations depend very strongly on the details of the input and fine-tuning 
of input and fitting parameters. 



Symmetry breaking 




FIG. 5: Three different results for the fan plot curves, with &i_n from set 1. Full line: from 



Tab. VII, dashed line: from Tab. VIII dotted line: from Tab. \X\ The bands shown here are the 



one-sigma error bands pertaining to the full line. Recall that the fits corresponding to the dashed 
and dotted lines include the experimental baryon masses. 



Let us now discuss the parameters which determine the SU(3) symmetry breaking vi- 
sualized by the fan plots, as seen in Fig. [5} First, it is reassuring; to observe that it does 
not matter whether bo,F are used as fixed input in the full fits or not: Tables 



VII 



and 



IX 



and Tables VIII and [X] are very similar, even though the input values in Tab. VII and |VIII" 
were obtained with approximations linear in the variable v for the ratios fs and a quadratic 
approximation for X N , and truncated chiral expansions for f B (eq. ([68]) ) , including a re- 
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ordering of terms to eliminate mo in favor of m* in the ratios. The bounds from the linear 
fits given in eqs. ( 75[76 ) are respected in all our resulting parameter sets. Even a rather 
large variation of the meson decay constant (from ~ 112 Me V to F ~ 80MeV) does not 



invalidate this picture, as can be seen from the fourth and the last column of Tab. XI It is 



also interesting to note that the combination — bp falls into the bound given in eq. (|37l) 



for bo,F for all fits from Tabs. VII X, and is slightly above that range for the fits with set 



1-3 from Tab. IXII From all these observations we can conclude that the determination of the 
LEC combinations parameterizing the linear part of the flavor symmetry breaking is very 
stable, thanks to the impressive new data leading to the fan plots. This is also illustrated 
by the error band and the three different fit curves shown in Fig. |5j 

It remains to discuss the parameters determined from the nonlinear symmetry breaking ef- 
fects. They appear at fourth chiral order in our formulae, which is the highest order we 
can compute exactly within the present one-loop framework. Therefore, we cannot expect 
a very high precision here. From our results, we extract the following allowed ranges for the 
remaining parameter combinations (in units of GeV -3 ): 



-0.05 <di < 0.15, 
0.00 <d 2 < 0.25, 
-0.10 <d A < 0.20, 
-0.75 <d 7 < -0.05 . 



(88) 
(89) 
(90) 
(91) 



We stress again that these estimates concern the numerical values of the renormalized LECs 



at n = 1150 MeV, and should be evolved by eqs. (52) when using infrared regularization at 
a different scale \i. 

Obviously, the purely linear fits lead to an underestimation of the uncertainty in d 7 : Here 
we had to shift the lower limit of our previous constraint, mostly due to the fit results shown 



in Tab. XI, compare eqs. (78) and (91). In Fig. 6 we show the constrained behaviour of 



the function X^[y\ which is mostly secured by the combination d-j. Here, the expectation 



of convergence at low orders must fail because the higher order terms in eq. (72) have to 



compensate for the leading (0(p 3 )) effect. This might partly explain why the range for dj 



in eq. (91) turns out to be comparatively broad. 



Including the experimental masses in the fit mainly (and strongly) affects the result for 
<i 2 , which parameterizes the curvature of the fan plot curves for the 5 and the nucleon, 
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compare e.g. Tab. VII and Tab. VIII , and the full versus the dashed and dotted curves in 
Fig. [5] While in all cases d% is of natural size, this shows that a reliable chiral extrapolation 
to the physical point, using one-loop BChPT expressions together with fan plot data, must 
still be considered problematic, in spite of the improvements in determining many LEC 
combinations worked out in the present contribution. In the framework of Heavy-Baryon 
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FIG. 6: The function X^{y) for the four fits of Tab. VII Red: set 1, green: set 2, blue: set 3, 



black: set 4. The band shown here is the one-sigma error band pertaining to the red line. 

ChPT, the low-energy constants d\-j have been estimated from the resonance saturation 
principle, combined with a fit to experimental masses and sigma terms, by Borasoy and 
Meifiner [I], see their eq. (68). Setting /3r = in those equations, and forming our preferred 
combinations, we see that their results for d\^j and all lie in the ranges given above (their 
G?6j however, is much smaller than in our tables, d^ M = -0.022 GeV" 3 ). Moreover, according 



to their estimates, b ^ M = 0.176 GeV -1 is somewhat above our allowed range (see eq. (75)), 
while their bf M = —0.344 GeV -1 again agrees very well with our typical fit results. One 
should note that such a comparison can only be of qualitative nature, for various reasons. 
Besides the differences of the schemes and the issue of scale dependence, it is not known 
how accurate the estimates from resonance saturation can probably be here, in particular 
because these LECs are not of "dynamical" nature (see also the discussion in sec. 5.2 of 
[XT]). What is more, the masses and 0^(0) = 45 MeV at the physical point have been input 
in the estimates of [I] , and we have already seen that in general the convergence properties 
of the chiral expansions are far from satisfying there due to the large m s -corrections. 
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Sigma term 



Rewriting our formula for the contact term contribution to the sigma term, eq. (48), in 
terms of our fitting parameters defined in eqs. fl74[81| ), we find 



(X6 
46' + 2b D + 2b F + 9QB md 6 + — {B Q 5m t ){3d x + 3d 2 - Ad 3 + 9d 5 + Qd 7 

(92) 



(93) 



This shows that we can write the sigma term resulting from our formulae as 

g 

o-ttTv(O) = a^ N (0)\ d3=d5=0 + -(2B 5m e )(2B m e )(9d 5 - Ad 3 ) 

« ^(0)| d3=d5=0 - \Ml{Xl - Ml){9d, - Ad,) , 

where the first term is written in terms of our fitting parameters. The unknown constants 
ds,ds parameterize the dependence on fh of the terms linear in Srri£. The latter dependence 



starts as soon as one includes chiral loops, i.e. at 0(p ), see e.g. eq. (68). Let us consider 



two examples to estimate the unknown parameter combination in question, ignoring the 



issue of convergence for the moment. For the fit of Table \ VTl\ with set 1, the first term is 
07rjv(O) \d s =d 5 =o — 67 MeV, while for the fit with set 4, it results in 47 MeV. For the remainder, 
we can only give a very rough estimate based on a consistency condition: Requiring that the 
terms ~ d 3 , d^ do not yield more than a 100% correction in the combinations bo,b' F gives 
bounds |c/ 3 | < 0.07 GeV~ 3 and \d 5 \ < 0.38 GeV~ 3 -0.41 < (d 5 - §4) GeV 3 < 0.41, which 
produces sigma terms between 39 ... 95 MeV for set 1 and 19 ... 75 MeV for set 4. 

However, it is to be noted that we do not see a clear sign of convergence for the sigma 
term at physical quark masses, so one should not overinterprete these last results. Basically 
we have to conclude that the sigma term could only be reliably determined from BChPT 
extrapolations when more data also for smaller average quark masses (so that e.g. M* ~ 
300 MeV) becomes available. 



VIII. CONCLUSIONS 



In summary, we have analyzed lattice data for octet baryon masses from (25] employing 
manifestly covariant BChPT with three dynamical flavors. We were able to give bounds on 
the numerical values of certain combinations of low-energy constants (LECs), which are, in 
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our opinion, more reliable and accurate than it was possible before data for fixed average 
quark mass became available in [2S]: As we have seen, the simulation strategy used in that 
work offers some advantages for this purpose. Eqs. fl75|76| ), ( [84p5p6] ) and d88p9p0|9l| 
should be considered as our main results. We hope that these bounds will be useful in 
the near future, e.g. when studying the chiral behaviour of other hadron observables with 
the same simulation strategy. Though we have fixed the running of the baryon masses on 
the two trajectories {rh = const., 5mi} and {fh,Sm£ = 0}, there are still two undetermined 
parameters, chosen here as d^ and d§, which would have to be fixed in order to make an 0(p 4 ) 
prediction for the sigma term. However, even if those parameters were fixed accurately, the 
uncertainty in the determination of the sigma term would still be very large due to the very 
slow convergence near the physical point. For an investigation of sigma terms, a two-flavor 
version of ChPT is certainly superior; we refer to [42J for a recent determination. In closing, 
we stress again that additional data points for lower average quark masses would be very 
helpful in order to reach a truly controlled chiral extrapolation. 
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